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1. Introduction 

It is a common folklore that the chiral phase transition takes place at high temperature and/or 
density in Quantum Chromodynamics (QCD). Meson masses are crucially influenced by the chiral 
condensate, and thus are very interesting observables in the chiral phase transition region. 

It is possible to extract hadron masses quantitatively in the lattice Monte-Carlo (MC) simula- 
tions at zero density. However, MC does not work well at high densities because of the notorious 
negative sign problem of the Dirac determinant. Therefore it is instructive to investigate high den- 
sity matter in the strong coupling lattice QCD (SC-LQCD). In SC-LQCD analyses with a mean 
field approximation (MFA), we can derive analytical expressions of the effective potential as a 
function of temperature {T) and chemical potential (/i), and hence we can avoid the sign problem. 
The numerical values of physical observables derived analytically in SC-LQCD should be repro- 
duced in lattice MC simulations at least in the strong coupling region. From this point of view, 
Refs. [jl]] provides an interesting comparison of strong coupling results in MC and analytic studies. 

Chiral phase transitions are governed by the effective potential, which also plays an essential 
role to obtain the hadron mass in SC-LQCD. The effective potential at finite temperature and den- 
sity has been extensively studied under various conditions For example, we have recently 
derived analytical expressions of the effective potential as a function of T and }JL in SC-LQCD, 
including the baryon propagating effects or l/g^ corrections 0. On the other hand, while hadron 
masses have been studied at zero temperature |§[ |8|, ^, the expression of meson masses as a func- 
tion of T and /i has never been derived before. 

In this proceedings, we derive an analytical expression of meson masses as functions of the 
chiral condensate, which is a function of T and /i in SC-LQCD for color SU(A'^c)- 



2. A brief summary of meson mass derivation 

We start from the action and partition function of lattice QCD with one species of staggered 
fermion(;^) in the strong coupling limit, where we omit the pure gauge plaquette terms (oc 1 /g^), 



/ exp 



(2.1) 

^/ = ^"L[e''XxU,{x)X,^(,-e-^{h.c)] , M,= , (2.2) 

^ X a=l 

where mo, Uq, Uj represent the current quark mass, the temporal and spatial link variables, respec- 



tively. We introduce the lattice chemical potential jj. following the procedure in Ref. []10|]. 

First we perform the path integral over spatial link variables Uj, and keep only the leading 
order terms in the l/d expansion. 



Z = / exp -S; + ^^±—^-moY^M, + ff{l/Vd) 



(2.3) 
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Next we introduce the chiral condensate a °^ {XX) through the bosonization procedure, 



Z ~ 



exp 



1 



^ OmxV^ ^ (xy ) 5mn Ony - "^ix^ " 



mn,xy 



Vm (xy) = 7^ L (5x+;,y + \-ly) ' "^nJ = m + a„x , 



AN, 



(2.4) 
(2.5) 



0=1 



where "m" or denotes the temporal lattice site, which takes an integer value in [\,N]. Now we 
can perform the path integral over staggered quarks, and obtain the following effective action of a 
up to a constant, 

I &a , 5eff[a] = ^ £ a^xV„ i(xy)5^„c7„y -£log/?x[a] . (2.6) 

mn,xy x 



Integral over ;^ and ?7o generates the interaction term log /?x[o'], where Rx[(y] is given as. 



/?x[al 



dUo{x) 



'l,x^c 











(2.7) 



where Ic represents an A/^^. x unit matrix, and ?7o is given in the temporal gauge as, 

t/o(x)=diag{^''«'W,---,^.'«"^W}. 



(2.8) 



Periodic and anti-periodic boundary conditions for gluons and fermions are respected in the present 
finite T treatment as found in the form of Eq. (^!7|), where we find different signs in the upper-right 
and the lower-left components. 

We decompose the chiral condensates a„(x) to the equilibrium value a and fluctuations 5a„(x), 
a„(x) = a + 5a„(x). The equilibrium value a is determined from the stationary condition of the 
effective action. 



_2Nc 1 dR^ia] 







(2.9) 



d R[o] 50n 

The meson mass is defined as the pole of the propagator for 5a, which is obtained from the second 
order variation of the effective action (5^5eff). The stationary condition Eq. (2^) ensures that the 
first derivative of /?x [o] is independent from the space-time point. Thus when we require the null 
average fluctuation condition £„5a„ = 0, the first derivative of Rx[o] in disappears. The 

17-1/ \S S ^^'nn\(y\ 

V„ (xy)d,„„-dxy- 



second order variation SiS'eff is found to be. 



mn,xy 



R[a] 



fir, fi-P) ^ 



5^Rx 



(2.10) 



(2) (2) 

It is necessary to evaluate Rmn to obtain meson masses. We derive Rmn at finite T and jj. in the 
following section, which has not been derived previously. 
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3. Evaluation of quark hopping in the temporal direction 



We utilize the formulation developed by Damgaard, Kawamoto, Shigemoto Faldt and 
Petersson and Nishida [Q]. First we reduce the {Nc xN) x (Nc x N) determinant in Eq. (2.7) to 
Nc X Nc determinant in the form 



dUo det 



(3.1) 



where the quark hopping kernel Xn is a functional of a. We assume that the number of temporal 
lattice sites N is even and its inverse N^^ is identified as the temperature T. Since is a function 
of XjM, it is enough to evaluate the Uo integral in equilibrium, and this has been done by utilizing the 
Vandermonde determinant technique combined with the recursion formula [^, Q]. It is also possible 
to perform the Uq integral explicitly by using the one link integral technique The obtained R^ 
reads, 

+ 2cosh- ^ 



R,{XN[a]=Y + Y-',ii 



(3.2) 



In this proceedings, the baryonic effects are considered for just a temporal direction, which is 
reflected on the second term in Eq (3.2). 

Next it is necessary to evaluate the quark hopping kernel X^v, which is found to be [||], 



Xn =B\^...N + B2---N-\ ■, B\...^N 











-e ^ 2m 



N.x 



(3.3) 



In equilibrium (5a„(x) 0, C7„(x) 
equilibrium value of Y and B as. 



a), recursion relation B„ = 2m(^^B„_i +B„-2 leads to the 



Y[o^o\ 

{cosh[(?i+l)£']/cosh£' (« = even) 
sinh [(« + 1 )E\ I coshS (n = odd) 



(3.4) 
(3.5) 



where E = sinh ' (mo + cj) denotes the one dimensional quark energy. By substituting Y in R^, we 
get the effective potential from the effective action (Eq. (p^), 

-^ir\h[{Nc+l)E{6)/T] 



=^eff(^)=5eff(a)/^: 



1 2N, 



2 d 



riog 



sinh[£'(a)/r] 



+ 2cosh 



(3.6) 



This has been obtained also in Ref. [Q] and included in the effective potential in Ref. [0]. The 
critical values Tc{lJL = 0), pLciT = 0) and the tri-critical point Ttcp can be extracted from the effective 
potential Eq. ( |3.6| ) and have been studied in Ref. iQ] in the chiral limit, 

(3.7) 
(3.8) 



r,(At=0) = J(A^, + l)(A^, + 2)/[6(A^, + 3)] , M.(r=0)~0.55 



rtcp = [J225N^ + 20d^{3N^ + 6Nc-4)-l5N,]/{20d) . 
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(2) K 

Finally, it is necessary to evaluate Rmn to obtain the meson mass spectrum(See Eq. ( |2.10D ). 
Here, the translational invariance simplifies the calculation. Since the system has translational 
invariance, and are cyclic invariant for the temporal indices. By using this translational 
invariance we obtain, 



dX, 



N 



d^X 



dOn 

'n 



dXK 



d^X 



2B. 



: 4B,/_iB/v-n'-i {n=N—\n — m\ 



(3.9) 
(3.10) 



It is found that Eqs. (|3J), (^) and ( [3.10| ) lead to. 



R 



(2) 



-4(B )2^'^ 



- - dR 

+ 4B„'-lB/V-,/-l^r— 



(3.11) 



in a straightforward calculation. In equilibrium, we can evaluate {dR/dXf^)\a in the following way. 



1 dR 



RdX, 



N 



1 dR [dX, 



RdOn 



dOn 



-1 



2N, 



1 



d 2Bn^i 



(3.12) 



where we used Eq. (2.9) and (3.4) in the final equality. Taking into account the even number of 
temporal lattice cite (N =even), and substituting Eq. (53) and (p.l2|) into (p.l 1[), we obtain. 



(2) 



mn 



2Nce'''"''cosh\{N-2m')E] , , , 

-a : — r-, 1 — (fn = \n — ml 

d cosh£sinh[£'/r] /? ^ ' ' 



(3.13) 



where does not include temporal indices, and here we do not bother to write its explicit form. 



4. Meson mass 



stituting Eq. ( |3.13| ) into ( [2. 10| ), we obtain. 



Now we can explicitly evaluate t he second variation of the effective action (Eq. (2.10)). Sub- 

_ 2Nc e'™'cosh [(A^ - 2m')E] 



mn,xy 



d coshE sinh[£'/r] 



5 any , (4.1) 



where we ignore the effect of A^/R by requiring the null average condition: l^„5a„ = 0. We 
perform the Fourier transformation of Eq. ( |4. 1[ ) by using the translational invariance, Y.n^i = 



yN m — > y ' for the second term in Eq. (14. lb, and we obtain. 



I5a«(k) 



£0,k 



2Nc/d {2Nc/d)o{o + mo) 



+ 



Y,jCoskj cos (0 + 2{o + moY + I 



(4.2) 



Here we introduce the prescription proposed in Ref. 

(co,k) = (/M,0) + (5v)7r, 



(4.3) 
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where (5v) is a "t/ + 1" dimensional vector which takes or 1 and originates from the tastes degrees 
of freedom. In this prescription, "M" is regarded as a meson mass. By putting the obtained inverse 
propagator for 5a£o(k) equal to zero, the meson masses are found to be, 

±coshM^(a;r,jU) = 2{o + mQ)(^-^o + +\ (4.4) 

d 

K = Y^cosdj G {-d,-d + 2,--- ,d-2,d} . (4.5) 

Here, we tentatively take the plus sign to obtain the real number meson mass, and regard taste 
effects K as meson species in the same way as in Ref. In this scheme with d = 3, "fc = 
—3, —1,1, 3" corresponds to "tt, p, a^" and "5" mesons, respectively. 

For a small current quark mass mo ~ 0, we obtain, My^ ~ l^om^, which may be regarded 
as the PCAC relation. When the temperature T and quark chemical potential \i are given, the 
chiral condensate o{T is obtained from the effective potential Eq. (^), then the meson mass 
Mk{o;T,ij.) is determined by Eq. (4^). In Fig. [|, we show T and jj. dependence of meson masses 
in the lattice unit. With a given T lower than the tri-critical point T^p (c.f. Eq. (3^)), the meson 
masses except for pion discontinuously drop when jj. passes through the critical value Hc, where the 
phase transition is first order (upper-left panel of Fig. [I]). At higher temperatures, Ticp <T<Tc, the 
phase transition is second order and thus the meson masses except for pion quickly but smoothly 
decrease to zero when jj. approaches to /Xc (upper-right panel of Fig. |l|). When T approaches to its 
critical value with a given /i lower than /i^, the meson masses quickly and smoothly decrease 
to zero, where the phase transition is second order (lower panels of Fig. |l]). The T dependence of 
meson masses is slightly affected by a given jj. in the chiral broken phase (lower panels of Fig. |l|). 



5. Summary 

We have investigated the meson mass spectrum at finite temperature and density by consider- 
ing the leading order of the l/d expansion in the strong coupling limit of lattice QCD. We have 
derived an analytical expression of meson masses as functions of the chiral condensate which is a 
function of temperature and chemical potential. We have thus explicitly studied the temperature 
and chemical potential dependence of meson masses near the critical value. We have confirmed 
that the pion mass satisfies the PCAC relation, and that masses of other mesons decrease to zero 
very quickly when the chemical potential or temperature approaches to the critical value. 
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